form  ^ppro^ccf 
QMS  A/o.  070^~C1B8 


lg/UU3 


REPORT  D0CUMEiMTAT50^J  PAGE 


-or  'fi-  >>:  ai-iu'c  ■ .  i  .u'  o-.t  ;hv 

_-r  t.'  ■•‘i.Mnl.tir'ir.  j  n,-_-Jcd.  j,r,j  rj  1  •,  on  •  .  J  ir. :  j'-'i  j L'V"  ■>  '■'iJ 

■'  — ■.:  ciV*.!  >'■ 'll!'’ n_  .r.  luLJ::'iy  jCT)*." ’■>*  -  e  -iiU'-  <>*  i  «'>-Mcu^nrr'.  jv^v-i.-.-v  r-r  .v^nf 

■  v  i.’n.;.  -'t.lf  *  •)'  0  t'?  .y  v n. .-<•  ..-rv-ni  .mu  fieOj-.L” 


0.'  J|..|-.r.',  .>.l»tin;3  C.lt.t  VlUfCTS 

t  -.ri-.  bu' •' -T  i»\T{'  yr  ^n,  v'ricr  ij!  xn'5 

i'r  ror  -dI  Cjo-.T-''<':  ■’iv  .nd  12  5?  3-^ifrr-.i^n 

.  r  ■'.* J.'ct  {u  ,»0^-0  '  yy*.  '.Vv-hir'giO'i  7-r  .-O VJ 1. 


■  AGePiCY  use  ONLY  (iejv/e  O/an/r)  h  I  DATE 


(  3  REPORT  type  AMO  DATES  COVERED 


i|  November  20,  1995  i _ Final.  30  Sep.  93  -  30  Sep.  9 


!■'  a.  title  and  s’jdTjTU 

r 

.1 

j'  Computational  Structures  of 
1  Multirate  Filter  Banks 

rimr^osisj  '  — — 

! 

j  Richard  Tolimieri 

A 


2?.  PEKEORMIKG  ORGAN52AT!ON  NAWE(5)  AMD  ADDRESS(SS) 

i  The  City  College  of  New  York 
Ij  Convent  Avenue  &;  138th  Street 
j|  New  York,  NY  10031 


I'o.  SPONSORING.  ViOWITORlNG  AGENCV  NAME(S)  AImO  aOOnESS(fS) 


5.  EUWOINC  NUMScRS 


F49620-93-1-0606 


■'  8.  PERFORMING  ORGANIZATION 
REPORT  number 

rBros'f^T^- 

j  C\Ls:^CC\'^ 


R  10.  SPOMSORIWG/MOWITORIWG 
!!  AGEMCY  report  NUtViBER 


I  AFOSR/NM 

|i  110  Duncan  Avenue,  Suite  B115 
Bolling  AFB  DC  20332-0001 


!  DISTRIBUTION /AVAILABILITY  STATEMENT 


12b.  DI5TRI8UTION  COOL 


:|  13.  ASSTRACT  {M^xirr^ urn  200  words) 

r  We  have  created  new  tools  for 

■  •  Developing  a  formalism  in  terms  of  shift  and  stride  permutation  operators  and  scalable 
'  tensor  product  kernels  for  modeling  complex  networks  of  multirate  filter  systems  and  an 
:  algebra  for  automatic  manipulation  of  fundamental  parameters  including  network  topologies. 

•  Applying  new  mathematical  tools  to  the  theory  of  multidimensional  multirate  filter 
structure  which  bring  out  the  fundamental  role  played  by  direct  sum  decompositions  and 
short  exact  sequences  and  permit  a  unified  treatment  of  both  standard  and  new  algorithms 
featuring  new  topologies,  communication  paths  and  data  structures. 

•  Developing  new  algorithms  for  constructing  families  of  unitary  transforms  (group  tr- 
nasforms)  from  the  Clifford  theory  of  idempotents  and  studying  the  applications  of  these 
unitary  transforms  to  multirate  filtering  and  transform  coding. 


Tl1%WS 


I  15.  NUMBER  Of  PAGES  g 


Multirate  Filters,  Tensor  Products,  Stride  Permutations  |  is  price  cqoe  — - 1 

Network  Topologies,  Short  Exact  Sequences,  Idempotents  i;  i 


-.7.  SECUS:'.'  .  CA;<;;?’C.AT;0N  s  ;G  security  CtASSlFlCATiON  r  is.  SSCURITV  classification  1;  20.  UMITATiOM  of  abstract  | 


il'.T.  SECU3..  - 
;■  Cr  FirAOT' 

Unclassified 

'i 


N<;m  75^U)-*J7-.:80  SSOO 


OF  ISIS  PAGC 

Unclassified 


OF  ABSTRACT 

Unclassified 


SrcnrqV.rH  F[5rm  298  2-89) 

HU 


19960207  023 


QTI.^-,!.rT¥  OTEFSCISE)  I, 


This  report  consists  of  three  parts:  research,  collaboration  and  academic. 


Research 

Rapid  advances  in  system  design  tools  and  computer  technologies  offer  the  poten¬ 
tial  of  parallel  and  distributed  computing  of  highly  complex  and  computationally 
intensive  multirate  filter  systems.  Typically  applications  to  image  processing,  com¬ 
munications  and  automatic  target  recognition  feature  the  need  to  process  and  analyze 
large  quantities  of  data  by  a  variety  of  sophisticated  mathematical  techniques.  Usu¬ 
ally  an  extensive  special  programming  effort  is  required  to  achieve  necessary  levels  of 
performance  which  must  often  be  repeated  for  even  small  changes  in  system  design 
or  target  architecture.  This  effort  can  be  extremely  time  consuming  and  error  prone 
greatly  limiting  and  delaying  the  impact  of  these  new  technologies  on  the  science 
and  engineering  community.  During  the  last  year  we  have  created  new  tools  which 
promise  to  utilize  and  extend  these  advances  at  minimal  time  and  cost  in  software 
development  by 

•  Developing  a  formalism  in  terms  of  shift  and  stride  permutation  operators  and 
scalable  tensor  product  kernels  for  modeling  complex  networks  of  multirate  filter 
systems  and  an  algebra  for  automatic  manipulation  of  fundamental  parameters 
including  network  topologies. 

•  Applying  new  mathematical  tools  to  the  theory  of  multidimensional  multirate 
filter  structure  which  bring  out  the  fundamental  role  played  by  direct  sum  de¬ 
compositions  and  short  exact  sequences  and  permit  a  unified  treatment  of  both 
standard  and  new  algorithms  featuring  new  topologies,  communication  paths 
and  data  structures. 

•  Developing  new  algorithms  for  constructing  families  of  unitary  transforms  (group 
trnasforms)  from  the  Clifford  theory  of  idempotents  and  studying  the  applica¬ 
tions  of  these  unitary  transforms  to  multirate  filtering  and  transform  coding. 

Some  aspects  of  this  research  will  be  presented  in  an  invited  address  at  102nd 
annual  meeting  of  the  Americal  Mathematical  Society  in  Orlando,  FL  in  January 
1996.  Moreover  this  research  has  resulted  in  several  papers  and  presentations  at 
conferences. 

•  ’’Algebra  of  Multidimensional  Multirate  Structures,”  submitted  for  presenta¬ 
tion. 

•  ’’Group  Algebras  and  Orthogonal  Decompositions,”  preprint. 
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•  ”Weyl-Heisenberg  Systems  and  Finite  Zak  transform,”  accepted  for  publication, 
Journal  of  Visual  Communication  and  Image  Representation,  Academic  Press, 
Inc. 

•  ’’The  Computation  of  Weyl-Heisenberg  Coefficients  for  Critically  Sampled  and 
Over-sampled  Signals,”  SPAT  ICSPAT94,  2  pp824-829,  Oct.  1994. 

I  expect  and  have  verbal  agreement  that  complete  results  of  this  research  will 
appear  in  bookform  to  be  published  by  Birkhauser. 

Several  student  projects  at  CCNY  have  implemented  codes  analyzing  and  testing 
the  efficiency  of  the  algorithms  developed  in  this  research  as  well  as  thier  applicability 
to  transient  signal  detection,  transform  coding  and  signal  analysis. 


Collaboration 

•  With  Professor  Patrick  Combettes  of  CCNY,  a  project  is  under  way  to  combine 
time-frequency  and  time-scale  analysis  methods  with  convex  projection  methods 
for  applications  to  large  scale  image  processing.  An  NSF  proposal  has  been 
submitted  with  professor  Robert  Alfano  of  CCNY  as  the  principal  investigator. 

•  With  Professor  A.J.  Devaney  of  Northeastern  University  a  project  is  underway 
to  develop  and  implement  wavelet  based  multirate  filter  algorithms  for  several 
projects  including  applications  to  ultreisound  tomography  and  radar  imaging. 

•  With  Professor  A.J.  Devaney  of  Northeastern  University  a  project  is  underway 
to  study  new  physical  models  and  algorithms  to  solve  phase  retrieval  problems  in 
imaging  and  crystallography.  One  aspect  of  this  research  in  collaboration  with 
Professor  B.  York  of  Northeastern  University  applys  to  the  problem  the  new 
single  isomorphic  replacement  methods  introduced  by  Professor  C.  Giacovazzo 
of  University  of  Bari,  Italy.  Several  papers  are  in  preparation  including  one 
which  develops  the  role  of  data  symmetry  in  reducing  and  organizing  code  for 
fundamental  triplet  calculations.  This  collaboration  had  resulted  in  submitted 
and  plans  for  future  proposals. 

•  Atlantic  Aerospace  Electronics  Corporation  has  been  developing  algorithmic 
strategies  for  the  efficient  execution  of  SPW^^  programs  and  other  system  de¬ 
sign  tools  in  a  multiprocessor  environment.  These  system  design  tools  have 
the  capability  of  creating  custom  coded  blocks  and  building  complicated  func¬ 
tionality  by  heirarchical  block-nesting,  but  tools  for  translating  complex  system 
designs  into  efficient  multiprocessor  codes  are  still  primitive.  With  Dr.  R.  Orr 
of  Atlantic  Aerospace  Electronics  Corporation  a  still  pending  SBIR  to  NSF 
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has  been  written  whose  goal  is  to  develop  sophisticated  tools  for  carrying  out 
efficient  multiprocessor  implementations. 


Grants  and  Contracts 

FY  95  AFOSR  MURI  Topic:  Computational  Electromagnetics 
Subcontract  Topic:  Mathematics  of  Computing 

STTR,  F49620-95-C-0073 

Thermal  Analysis  of  Multichip  Modules  Using  Domain 
Decomposition  and  Wavelet-Capacitance  Matrix 


Academic 

The  graduate  program  in  Electrical  Engineering  at  CCNY  has  been  subject  to  much 
criticism  during  the  last  few  yeaxs  as  tuition  increases  and  the  difficulty  of  post  gradu¬ 
ate  employment  has  greatly  reduced  the  student  body  with  a  subsequent  reduction  in 
cost  offerings.  The  program  has  a  large  minority  enrollment  having  limited  network¬ 
ing  skills  as  to  relevant  fields  of  study  and  as  to  possible  employment  opportunities. 
One  purpose  of  this  sabbatical  was  to  bring  to  these  students  a  more  focused  program 
of  digital  signal  processing  based  on  an  understanding  of  those  fundamental  courses, 
engineering  culture  and  application  areas  which  would  best  serve  their  intellectual  and 
career  potential.  A  related  goal  was  to  establish  contact  and  possibly  close  relation¬ 
ships  with  potential  sources  of  problems  and  employment  so  that  students’  research 
efforts  could  be  better  directed.  Multirate  filtering  and  multiprocessor  algorithmic 
development  were  chosen  as  target  areas  as  these  tools  are  expected  to  be  of  long  time 
significance  to  a  wide  range  of  applications  to  the  DoD  and  medical  imaging  com¬ 
munity.  During  this  last  year  important  contacts  have  been  made  at  Northeastern 
University,  Hanscom  Air  Force  base  and  Atlantic  Aerospace  Electronics  Corporation 
which  have  already  born  fruit  as  to  employment  opportunities.  Moreover,  at  CCNY 
several  new  courses  in  DSP,  multirate  filter  design  and  computer  architecture  have 
been  introduced  and  serve  as  a  basis  to  new  courses  and  proposals  in  imaging,  radar 
and  multiprocessor  implementation. 


Publications  during  the  Grant  Period 

1.  ’’Algebra  of  Multidimensional  Multirate  Structures,”  submitted  to  Journal  of 
Applied  Mathematics,  North-Holland. 
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2.  ’’Group  Algebras  and  Orthogonal  Decompositions,”  preprint. 

3.  ”Weyl- Heisenberg  Systems  and  Finite  Zak  transform,”  accepted  for  publication, 
Journal  of  Visual  Communication  and  Image  Representation,  Academic  Press, 
Inc. 

4.  ’’The  Computation  of  Weyl-Heisenberg  Coefficients  for  Critically  Sampled  and 
Over-sampled  Signals,”  SPAT  ICSPAT94,  2  pp824-829,  Oct.  1994. 

5.  with  R.  Orr,  “Poisson  Summation,  the  Ambiguity  Function,  and  the  Theory  of 
Weyl-Heisenberg  Frames,”  The  Journal  of  Fourier  Analysis  and  Applications  1 
3,  1995. 

6.  with  A.  Brodzik,  “Weyl-Heisenberg  Systems  and  the  Finite  Zak  Transform,” 
accepted  for  publication. 

7.  with  I.  Gertner,  “Multiplicative  Zak  Transform”,  Jour.  Visual  Comm,  and  Im. 
Rep.  6  1,  March  1995,  89-95. 

8.  with  M.  An,  Y  Abdelatiff  and  N.  Anupindi,  “Group  Invariant  Fourier  trans¬ 
form,”  Adv.  in  Electronics  and  Electron  Phys.,  93,  P.  Hawkes,  Ed.,  Academic 
Press,  1995. 

9.  With  G.  Kechriotis,  M.  An,  M.  Bletsas  and  E.  Manolakos,  “A  New  Approach 
for  Computing  Multidimensional  DFT’s  on  Parallel  Machines  and  its  imple¬ 
mentation  on  the  iPSC/860  Hypercube,”  IEEE  Trans,  on  Signal  Proc.  431, 
272-285,  January,  1995. 

10.  with  S.  Qian,  C.  Lu  and  M.  An,  “Self-Sorting  In-Place  FFT  Algorithm  with 
Minimal  Working  Space,”  IEEE  Trans,  on  Signal  Proc.  4210,  2835-2836,  Oct- 
ber,  1994. 

11.  with  M.  An  and  C.  Lu,  “DSP  algorithm  design  and  implementation  on  RISC 
architectures,”  presented  at  the  first  international  conference  on  Electronics  and 
Information  Technology  (ICEIT’94)  held  in  Beijing,  China. 

12.  with  M.  An  C.  Lu  and  S.  Qian,  “Self-sorting  In-place  FFT  Algorithm  with 
Minimum  Working  Space”,  IEEE  Trans,  on  Signal  Processing,  4210,  1994. 

13.  with  M.  An  and  C.  Lu,  Mathematics  of  Multidimensional  Fourier  Transform 
Algorithms,  October  1993,  Springer- Verlag,  New  York. 

14.  with  M.  An,  C.  Lu  adn  S.  Qian,  “A  Hybrid  Parallel  FFT  Algorithm  without 
Interprocessor  Communication,”  Proc.  of  IEEE  Inter.  Conf.  on  ASSP,  1993. 
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